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We consider dynamics of spatial beams in a dual-channel waveguide with competing cubic and 
quintic (CQ) nonlinearities. Gradually increasing the power in the input channel, we identify four 
different regimes of the pulse's coupling into the cross channel, which alternate three times between 
full pass and full stop, thus suggesting three realizations of switching between the channels. As 
in the case of the Kerr (solely cubic) nonlinearity, the first two regimes are the linear one, and 
one dominated by the self-focusing nonlinearity, with the beam which, respectively, periodically 
couples between the channels, or stays in the input channel. Further increase of the power reveals 
two novel transmission regimes, one characterized by balance between the competing nonlinearities, 
which again allows full coupling between the channels, and a final regime dominated by the self- 
defocusing quintic nonlinearity. In the latter case, the situation resembles that known for a self- 
repulsive Bose-Einstein condensate trapped in a double-well potential, which is characterized by 
strong symmetry breaking; accordingly, the beam again abides in the input channel, contrary to an 
intuitive expectation that the self-defocusing nonlinearity would push it into the cross channel. The 
numerical results are qualitatively explained by a simple analytical model based on the variational 
approximation . 
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^ '. I. INTRODUCTION AND THE MODEL 

^0 ' Dynamics of spatial light beams in nonlinear planar waveguides with a multi-channel built-in structure is a subject 
Ph , of fundamental interest to nonlinear optics, and has vast potential for applications to all-optical routing of data 
. streams. The basic model of the multi-channel system amounts to the nonlinear Schrodinger (NLS) equation for the 
• ^ ' local amplitude '^'{z^x) of the electromagnetic wave, where z and x are the propagation and transverse coordinates 
in the waveguide. In a normalized form, the NLS equation is |3| 
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j^-. + ^.. -t/(x)* + (5n(|*|2)* = 0, (1) 



where the second term accounts for the diffraction of light in the paraxial approximation, U{x) is the profile of the 
transverse modulation of the refractive index that accounts for the multi-channel structure, and ^n(|^'p) is a nonlinear 
' correction to the refractive index. In the case of the ordinary self-focusing (SF) Kerr nonlinearity, Sn = n2|^'P, 

lO with 7^2 > 0. 

The interplay of the diffraction and SF nonlinearity gives rise to channel-trapped spatial solitons. In the model with 
the Kerr nonlinear term and harmonic transverse modulation, U{x) = Uq cos (27ra;/L), trapped solitons and a possibil- 
ity of their switching were studied in Ref. It was found that the soliton's integral power, P = |^'(a;)|2(ix,which 
is a dynamical invariant of the NLS equation, may take any value, < P < oo, the entire solution family being stable. 
Switching, i.e., transfer of the soliton into an adjacent channel, was demonstrated under the action of a "hot spot" 
modeled by an extra term ~ S{x — xo)S{z — zq)^> in the equation, with xq fixed at the midpoint between the channels. 
J> The same model as in Ref. 0|, with z replaced by time t (but without the hot-spot term) was later introduced 

for an effectively one-dimensional Bose-Einstein condensate (BEC) with attraction between atoms, loaded in an 
optical-lattice potential j^. Spatial solitons in a model of the Kerr-nonlinear optical waveguide with the periodic 
' transverse modulation in the form of the Kronig- Penney profile (a periodic array of rectangular troughs, see below) 
were considered in Refs. Q. In optics, a system of parallel guiding cores can be created not only through a small 
periodic variation of the material's refractive index, but also in a virtual form, by means of a system of strong 
transverse laser beams illuminating a sample of a photorefractive material in the ordinary polarization, while the 
signal beam is launched in the extraordinary polarization . 

The functionality of optical beam-handling schemes may be greatly enhanced in media with competing nonlinear- 
ities, the simplest example of which is provided by a combination of the SF cubic and self-defocusing (SDF) quintic 
terms. Observations of the cubic-quintic (CQ) optical nonlinearity were reported in the PTS crystal 5], chalcogenide 
glasses and some other materials 0. In fact, the CQ nonlinearity occurs due to an intrinsic nonlinear resonance 
in the material, which also gives rise to strong two-photon absorption however, experiments with spatial solitons 
require short propagation distances, on the order of a few centimeters ^ |^, over which effects of the loss may be 
insignificant |10l| . 

The NLS equation with the CQ nonlinearity and the multi-channel guiding structure can be rescaled into the 
following form ,1IJ: 

+ = C/(a;)*-2|*|2* + 1*1%. (2) 
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In the uniform medium [U{x) = 0], Eq. ^ gives rise to a well-known family of stable solitons ^ 

^'(x, z) = exp (ikz) 
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(3) 



with the propagation constant k taking values in a finite interval, < /c < 3/4. In Ref. [ill, solitons were studied 
m CQ model (0) with a single guiding channel of a rectangular shape. A remarkable feature of the channel-trapped 
CQ solitons is bistability of the guided-soliton family: in the same region as mentioned above, < A: < 3/4, the 
channel supports a single soliton state, but the full existence region extends to an additional interval, 3/4 < fc < fcmax, 
where two different solitons are found for each k, both being stable (the soliton bistability does not occur in the same 
guiding channel if the model contains only the cubic nonlinearity [13 )• The CQ nonlinearity was combined with the 
multi-channel structure of the Kronig-Penney type in Ref. |14| . In addition to fundamental solitons, many families 
of higher-order solitons, also featuring the bistability, were found in the latter model. 

For applications to all-optical switching, more relevant solutions are not static solitons, but rather dynamical 
states featuring periodic transfer (coupling) of the beam between two adjacent guiding cores. In this context, the 
most appropriate model is one with two channels, similar to nonlinear directional couplers operating in the temporal 
domain, which are based on dual-core fibers (see Refs. and references therein). The objective of the present work 
is to report numerical results and some analytical approximations for the beam switching in the dual-channel CQ 
model, which is based on Eq. ^ with 

_ r 0, < iL and >^L + D, 

where D, Uq and L are, respectively, the width and depth of each channel, and the thickness of the buffer layer 
between them. A fundamental difference from earlier studied two-channel models with the Kerr nonlinearity is 
expected dependence of the coupling on the beam's input power P: in the Kerr model, the increase of P leads 
simply to transition from the nearly-linear-coupling regime to a nonlinear one, in which the beam stays in the input 
channel (coupling suppression). In the CQ model, a more sophisticated dependence on P may be anticipated, due 
to the competition between the cubic SF and quintic SDF terms. The main qualitative result demonstrated below 
is that the coupling suppression is changed by restoration of the coupling at larger values of P (at which the SF 
and SDF terms are in balance), and finally again by suppression of the coupling at very large P, when the quintic 
term dominates, even if the suppression of coupling by the seli-defocusing nonlinearity is a counterintuitive effect. In 
Section m we display typical examples of the four distinct transmission modes, and then present a diagram which 
provides for overall description of the switching in the dual-channel waveguide with the CQ nonlinearity. Using values 
of strengths of the cubic and quintic SF and SDF nonlinearities in available materials, we estimate that the switching 
can be achieved at power densities of the input beam on the order of GW/cm^, the corresponding coupling distance 
being ^ 1 cm. 

In addition to the numerical results, in Section IIIII we report results of a simple analytical approximation that 
provide for a qualitative explanation of the switchings in the present model. The analytical approach is based on a 
variational approximation, similar to that developed earlier in detail for temporal signals (solitons) in dual-core fibers 
[11 [HO (see also Chapter 6 of review P|). 

II. COUPLING DYNAMICS AND SWITCHING DIAGRAMS 

To simulate the transmission of the spatial beam in the system described by Eq. (jSJ, we integrated this equation 
with initial conditions (at z — 0) corresponding to a Gaussian beam, with amplitude A and width W, launched into 
one channel: 



*o = ^ exp 



VU2 
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(5) 



[x = {L + D) /2 \s the center of the first channel as per Eq. Q]. The splitting of light between the two channels is 
quantified by respective integral powers, 

(L + D)/2 -L/2 

Pi= j P2= j 1*1' da;. (6) 

L/2 -{L+D)/2 
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Equation ^ was solved numerically by means of the split-step method. The input power was gradually increased, and 
the evolution of the power distribution across the waveguide was monitored in the course of the simulations. Collecting 
numerical results, it was found that they are adequately represented by fixing the parameters at L = D = Uo = W — 2. 
This implies that the width of the channels is equal to the thickness of the buffer layer between them, and the diffraction 
and channel-trapping terms are comparable to the nonlinear ones when they are in balance. Results for other values 
of the parameters will be displayed too. 

If the input power is small, nonlinear terms in Eq. Q are negligible in comparison to the diffraction, which leads 
to the ordinary picture of the periodic coupling of power between the channels, as shown in Fig. 1. Note that the 
transmission length equal to coupling length .^coupi, i-e., half-period of the coupling (as in Fig. 1), provides for the 
full transfer of the beam into the cross channel. 

Increase of the input power makes the SF cubic nonlincarity dominant, switching (as usual) the transmission mode 
from the periodic coupling to one with the power abiding in the input channel. This situation is illustration by Fig. 
2. 

The two transmission modes shown in Figs. 1 and 2, and the switching between them are qualitatively the same 
as in ordinary dual-core models with the Kerr nonlincarity 151- Switching to an essentially new transmission regime, 
which is specific to the CQ medium, occurs with further increase of the power, in a range where the SF cubic and SDF 
quintic terms balance each other, rendering the situation similar to that in the linear model. This new quasi-linear 
regime is illustrated by Fig. 3. It is noteworthy that the linear-coupling mode in Fig. 1, and its strongly nonlinear 
but similar counterpart in Fig. 3 have nearly equal coupling lengths: for equal values of the system's parameters and 
beam's width, zt-f-V^^S « 7/6. 

Extra (final) switching occurs, if one keeps increasing the input power, to a transmission regime in which the SDF 
quintic term dominates. "Naively" , one might expect that SDF nonlincarity would only enhance coupling between the 
guiding cores; however, it actually suppresses the coupling, as shown in a typical example in Fig. 4. As observed in 
the figure, this regime is not completely equivalent to one in which the coupling was suppressed by the dominant SF 
cubic term (cf. Fig. 2), as the present transmission mode features conspicuous oscillations; nevertheless, the coupling 
is largely suppressed in this case too. 

A result similar to that demonstrated in Fig. 4 is known in a model of a Bose- Einstein condensate (BEC) with 
repulsion between atoms (which is formally tantamount to SDF cubic nonlincarity) trapped in a double-well potential, 
the two potential minima playing the same role as the dual-channel configuration. The corresponding Gross-Pitaevskii 
equation (the counterpart of the NLS equation for the BEC) gives rise to stable asymmetric states (see Ref0| and 
references therein). This symmetry-breaking effect has also been observed in experiment 0| [in Ref. a similar 
effect was demonstrated, theoretically and experimentally (using a photorefractive crystal), in a medium with saturable 
nonlincarity] . 

Results of systematic simulations are summarized in the form of a switching diagram, which is displayed in Fig. 5. 
To measure the coupling rate, we define s = (^2)inax /-^' maximum share of the total input power (P), which 

is coupled into the second channel. The diagram plots e versus the total power itself; to demonstrate the generic 
character of the dependence, it is presented for three different values of the width of the input pulse © (including 
W = 2, for which examples of the transmission modes are displayed in Figs. 1-4). The regimes with e approaching 
1, and with e <C 1 are, respectively, those featuring nearly complete periodic transfer of the power between the 
channels, and suppression of the coupling. In accordance with what was displayed above by dint of examples, the 
strong coupling takes place at small P, which corresponds to the nearly linear transmission, and in a relatively narrow 
interval of P providing for the balance between the SF and SDF nonlinearities. In two broad intervals, where either 
the cubic term or its quintic competitor dominate, the coupling is suppressed. For applications, a beneficial feature of 
the £{P) dependences in Fig. 5 is the steep transition between regimes of the two types, i.e., sharp switching. Figure 
5 also suggests that exact values of the power at which the switchings take place can be tuned by adjusting the width 
of the input pulse. 

To conclude this section, we note that, with typical values of the cubic and quintic nonlinear coefficients for 
chalcogenide glasses reported in Refs. ^2 = 2.2 x 10^^ cm^/GW, and = -7.9 x 10^^ cm'*/(GW)^ (they are 
defined so that the total nonlinear correction to the refractive index is 6n = + n^I^ , where / is the power density 
in the optical beam), the switchings specific to the CQ medium are attained in a ballpark of / = 3 GW/cm^, which 
provides for the balance between the SF and SDF effects, and is achievable in the experiment Q. Analysis of the 
corresponding wave-propagation equation written in physical units demonstrates that the respective coupling length 
(in the quasi-linear regime of balanced SF and SDF) may take values in the range of 0.1 to 1 cm. 
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III. ANALYTICAL APPROXIMATION 



The underlying equation (0) can be derived from the Lagrangian 



L = 



dx. 



(7) 



To apply the variational approximation to the present model, we assume a situation with very narrow and deep 
channels, separated by a broad buffer layer. The field configuration is then approximated as follows: inside each 
channel, i.e., in regions \x- {L± D) /2| < D/2 [see Eq. (g))]. 



'^{x,z) = yl±(z)cos 



x-{L±D)/2 
D 



(8) 



with complex amplitudes A± (this ansatz emulates the ground state of a quantum-mechanical particle in an infinitely 
deep potential box, therefore it vanishes at edges of the channel) . The form of the ansatz outside the channels is also 
borrowed from quantum mechanics, emulating a superposition of exponentially decaying ground-state wave functions 
behind the edges of very deep but finite potential boxes. 



^{x,z) = ^ A±(z)exp (~^Uo-{^lD) 



L±D 



(9) 



The substitution of the inner and outer parts of the ansatz, © and (Q, in Lagrangian Q and integration lead to the 
following effective Lagrangian: 
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{L + D)] {A+A*_+A^Al) , 



(10) 



where the prime stands for d/dz. Note that the last term in H10|l . which accounts for the interaction between the 
beams trapped in the two channels, was calculated by means of a general method elaborated in Ref. (see also 
section 2.3 in review [3). This method only assumes that, in a vicinity of each channel, the tail created by the other 
beam is much weaker than the field of the beam trapped in the given channel (note that L + D is the distance between 
centers of the two channels). 

The Euler-Lagrange equations, dLcs/dAj_ — dLcs/d [A'^Y 

R±e^'^±, where R± and are real amplitudes and phases: 

i?^ = kR^ sin (j), R'_ 



— can be cast in a real form, by defining A± 
sin 0, 



[R^Y + K {cos <j>):^, 



(11) 
(12) 



where the effective constant of the tunnel coupling between the two guiding cores is 



D 



\DJ 



\D 



(13) 



An obvious consequence of Eqs. (|ll|l and H12I) is the conservation of P = R\ + R^_ (it is proportional to the total 
power of the beam), and the fact that an equation for 0^ + 0_ is detached from the other equations. Thus, defining 



i?4 



Pcos{9/2),R- = VPsin(6i/2) 



(14) 



A = (P/8) (12 - 5P) , 



(15) 
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we end up with a dynamical system 



- = -Ksm0, 
(f) — X cos 6 — K (cos (f)) cot 61, 



which conserves the Hamiltonian, 



77= ^KCOS(/) — -Asin6'^ sin0. (17) 

Dynamical system p6|l is similar to that which was derived, in the simplest version of the variational approximation, 
for a two-component temporal soliton in a dual-core optical fiber with the Kerr nonlinearity [l^ Il7j| . The only 
significant difference is that the present CQ model gives rise to nonlinearity coefficient A which changes its sign with 
the increase of the beam's power P, see Eq. IjlSf) [k is always positive, according to its definition (|13|l ]. 

In the framework of this approximation, solutions corresponding to various coupling modes found above by means 
of numerical simulations of Eq. ^ with initial condition (O are obtained from the initial condition with the beam 
launched into one channel, i.e., 6' = or = tt, according to Eqs. 1)141) . As seen from Eq. (|17|l . this initial condition 
yields H = 0, hence the dynamical trajectory in the {9, (j>) plane, H = const, is generated by the equation 

cos(/) — ainO. (18) 

It follows from Eq. (|18|1 that the trajectories take completely different forms in cases |A| < 2k and |A| > 2k, 
see examples shown in Figs. 6(a) and (b). The former case corresponds to weak effective nonlinearity, i.e., to the 
coupling mode with either small input power, or with balanced SF and SDF terms, cf. Figs. 1 and 3, respectively. In 
qualitative agreement with the results of numerical simulations of Eq. in this (quasi- linear) case the trajectories 
run along the 6 direction, which means periodic transfer of power between the channels. Note that Eq. (|15() makes 
it possible to predict intervals of power P in which the quasi-linear coupling occurs: P < (2/5) (3 — V9^-"20k), and 
(2/5) (3 + y/9 — 20k) < P < (2/5) (3 + y/9 + 20k) [these expressions imply constraint k < 9/20, which comphes with 
the fact that expression (|13|l yields small k]. On the other hand, in the case of strong nonlinearity, |A| > 2k, i.e., 
in intervals (2/5) (3 - V9^^20^ < P < (2/5) (3 + V9 - 20k) and P > (2/5) (3 + V9-(-20k), Fig. 6(b) displays 
trajectories confined to a vicinity of = or 6* = tt, which corresponds to the suppression of coupling, cf. Figs. 2 
and 4. Thus, the variational approximation provides for a qualitatively correct explanation of the numerical results 
presented above. 



IV. CONCLUSION 



We have investigated possibilities of realizing power-controlled switching of spatial light beams in the dual-channel 
planar waveguide with the cubic-quintic nonlinearity. The interplay of the transverse diffraction and self-focusing 
(cubic) and defocusing (quintic) nonlinearities gives rise to two transmission modes similar to those known in models 
with the cubic nonlinearity, i.e., periodic coupling in the nearly linear situation, and suppression of the coupling in the 
mode dominated by the cubic term, and two additional regimes, one with the periodic coupling restored due to the 
balance between the competing nonlinearities, and the other (final) transmission mode with the coupling suppressed 
again when the self-defocusing quintic term dominates. Three switchings (between the four distinct transmission 
regimes) are not strongly affected by variation of the system's parameters, but they feature sharpness with variation 
of the power, especially the reverse switching from the coupling-suppressed mode dominated by the cubic nonlinearity 
to the strong-coupling one provided by the balance between the self- focusing and defocusing nonlinearities. A simple 
analytical approximation, which reduces the model to a Hamiltonian dynamical system with one degree of freedom, 
qualitatively explains all these coupling models and switchings between them. 
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No. 7001705), and by the Israel Science Foundation through a Center-of-Excellence grant No. 8006/03. R.D. and 
B.A.M. appreciate hospitality of the Optoelectronics Research Centre at the Department of Electronics Engineering, 
City University of Hong Kong. 
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Figure Captions 

Fig. 1. (Color online) The first half-period of the quasi- linear coupling regime, for input beam Q with A = 0.2 
and W = 2: (a) contour plots of the power distribution; (b) integral powers in the channels, defined as per Eq. (|^, 
vs. the propagation distance. 

Fig. 2. (Color online) The same as in Fig. 1, but after the switch to the transmission regime dominated by the 
self- focusing cubic term. In this case, the input amplitude is A = 1, i.e., the total power is larger by a factor of 25 
than in the linear regime displayed in Fig. 1. A longer transmission distance than in Fig. 1 is shown here, to make 
it sure that the coupling to the the cross channel never commences. 

Fig. 3. (Color online) The same as in Fig. 2, but after the reverse switch to the regime in which the self-focusing and 
defocusing terms (cubic and quintic ones) are in balance, thus again allowing strong coupling between the channels. 
In this case, the input amplitude in Eq. (jsj is yl = \/2, i-c, the total power is twice as large as in the case shown in 
Fig. 2. Unlike Fig. 1, about 1.5 periods of the coupling is shown here, to demonstrate the stability of the periodic 
coupling (actually, it remains stable indefinitely long, in more extensive simulations). 

Fig. 4. (Color online) The same as Fig. 3, but after the final switch to the transmission mode with suppressed 
coupling. The amplitude of the input pulse is ^ = 1.8, hence the total input power exceeds that in Fig. 3 by a factor 
of 1.62. 

Fig. 5. (Color online) Switching diagrams, which show the maximum share of the input power coupled into the 
cross channel, as a function of the input power itself. The diagrams are presented for initial Gaussian beams lO with 
three different values of the width, W — 1, 2, and 2\/2 (recall examples of the four transmission regimes displayed in 
Figs. 1-4 pertain to = 2). 

Fig. 6. Typical dynamical trajectories in the phase plane {9, (f)), as generated by Eq. H18|) : (a) A/(2k) = 0.5, which 
corresponds to the quasi-linear couphng mode; (b) A/(2k) = 2, corresponding to the coupling-suppressed mode, 
dominated by the nonlinearity. 
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